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We show how mapping techniques inherent to N 2 -dimensional discrete phase spaces can 
be used to treat a wide family of spin systems which exhibits squeezing and entanglement 
effects. This algebraic framework is then applied to the modified Lipkin-Meshkov-Glick 
(LMG) model in order to obtain the time evolution of certain special parameters related 
to the Robertson-Schrodinger (RS) uncertainty principle and some particular proposals 
of entanglement measure based on collective angular-momentum generators. Our results 
reinforce the connection between both the squeezing and entanglement effects, as well as 
1 allow to investigate the basic role of spin correlations through the discrete representatives 

of quasiprobability distribution functions. Entropy functionals are also discussed in this 
context. The main sequence correlations i-> entanglement i— > squeezing of quantum effects 
embraces a new set of insights and interpretations in this framework, which represents 
an effective gain for future researches in different spin systems. 

Keywords: Spin squeezing; entanglement; finite-dimensional discrete phase spaces. 

1. Introduction 

Originally introduced by E. Schrddinger in his seminal work on probability relations 
between separated systems^ "entanglement" indeed corresponds to a fundamental 
concept in phy sics that lies at the heart of many conceptual problems in quantum 
mechanics!^'' The mere comprehension of this abstract concept and its respective 
changing of status, in the recent past, for experimental measure, certainly represents 
a concatenation of efforts with significant progress in theoretical and experimental 
physics. It is important to stress that 'quantum entanglement' plays an essential 
role in multipartite quantum systems SI since its underlying physical properties can 
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be used as a specific resource for determined quantum-information taskP - re- 
cently, different proposals of quantumness correlations have appeared in current 
literature,^ emphasizing the possible connections with quantum entanglement.^ 
Therefore, understanding, identifying, measuring and, consequently, exploiting gen- 
uine quantum effects (or "quantumness"^) in multipartite systems constitute a set 
of obligatory prerequisites for grasping the fundamental implications of any quan- 
tum theory!^ ___ 

Recent theoretical proposali^ and experiment J ^ ^ I ^ involving measurements 
upon collective angular-momentum generators in different physical systems fulfil, 
in part, the aforementioned prerequisites, as well as corroborate the subtle match 
between spin-squeezing^^ and entanglement effects. However, some necessary ques- 
tions concerning the correlations among different spin components as chief agents 
responsible for spin-squeezing effects deserve be properly answered. Then, it seems 
reasonable to adopt a particular algebraic framework that, within all the inherent 
mathematical virtues, allows to: (i) map the kinematical and dynamical contents of 
a given spin system with a finite space of states into a N 2 -dimensional discrete phase 
space; (ii) comprehend the role of correlations through discrete Wigner functions 
and their connections with spin-squeezing effects; and finally, (iii) gain new insights 
on the effects exhibited in the Venn diagram below which lead us to take advantage 
operationally of their quantumness in many-particle experiments. Notwithstanding 
the appreciable number of papers in current literature propo sing similar theoretical 
frameworks with different intrinsic mathematical properties)^ let us focus our at- 
tention upon the formalism developed in Refs. 15-24 which complies such demands. 




Based on the technique of constructing unitary operator bases initially for- 
mulated by Schwingerf^l this particular discrete quantum phase-space approach 
embraces a well-defined algebraic structure where the N- dimensional pre-Hilbcrt 
spaces equipped with a Hilbert-Schmidt inner product!^ endorse the finite space 
of states. Moreover, it leads us (i) to exhibit and handle the pair of complemen- 
tary variables related to a specific degree of freedom we are dealing with (as well 
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as to recognize the quantum correlations between them), and also (ii) to obtain 
additional quantum information about the physical system through the analysis of 
the corresponding discrete quasiprobability distribution functions. Henceforth, the 
basic idea consists in exploring this mathematical tool in order to study the role of 
those correlations in connection with spin-squeezing and entanglement processes. 

Initially focused on some mathematically appealing features inherent to the col- 
lective angular-momentum generators, the first part of this paper basically discuss 
the spin coherent states and their relations with unitary transformations through 
a constructive point of view. Next, we establish a mod(N)-invariant operator basis 
which can be immediately employed in the mapping of quantum operators (acting 
into that particular N-dimensional space of states) onto well-behaved functions of 
discrete variables by means of a trace operation!^ These functions correspond to 
the representatives of the operators in a N 2 -dimensional phase space labeled by a 
pair of discrete variables for each degree of freedom of the physical system under 
investigation. Thus, all the necessary quantities for describing its kinematical and 
dynamical contents can now be promptly mapped one-to-one on such phase spaces. 
For completeness sake, we also present two different prescriptions that essentially 
determine the time evolution for both the discrete Wigner and Weyl functions. 

These results are then applied to an extended version of the LMG modelf^ 
which was initially introduced by Vidal and coworkers with the aim of investi- 
gating the statistical-mechanical properties and entanglement effects of a particular 
interacting spin system. Summarizing, the modified LMG model describes a finite 
set of spins | with a mutual anisotropic (XY) ferromagnetic interaction, and also 
subjected to a transverse magnetic field. The second part of this paper shows how 
the quantum correlations (here controlled by parameters related to the transverse 
magnetic field and anisotropy) affect the connection between the spin-squeezing 
and entanglement effects. To develop this specific task, we first employ the RS 
uncertainty principle^! in order to determine a spin-squeezing measure which in- 
corporates, in its definition, the covariance function — being such a function respon- 
sible for introducing, within the Heisenberg uncertainty principle, important contri- 
butions associated with the anticommutation relations of the angular-momentum 
generators. The subsequent comparison with the entanglement measure proposed 
in Rcf. [12] establishes, in this way, the aforementioned link and reinforces the fun- 
damental role of correlations in our initial analysis. Moreover, we also study the 
behaviour of the discrete Wigner function for certain particular values of time where 
occur an 'almost perfect match' between both the spin-squeezing and entanglement 
measurements (such procedure leads us to comprehend the intricate role of corre- 
lations in a finite-dimensional phas e sp ace). We have finalized this work with the 
analysis of some entropy functionals!^ For completeness sake, it is worth stressing 
that we have worked with a number of spins which preserves basic quantum effects, 
contrary to the thermodynamic limit. 

This paper is structured as follows. In Section 2, we present an important set of 
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essential mathematical tools related to the angular-momentum generators and spin 
coherent states, as well as an interesting discussion on unitary transformations and 
their corresponding geometric interpretations. In Section 3, we introduce a map- 
ping kernel for finite-dimensional discrete phase spaces whose inherent properties 
allow, within other important virtues, to describe the kinematical and dynamical 
contents of a given physical system with a finite space of states. In Section 4, we 
apply the quantum-algebraic approach developed in the previous section to the 
modified LMG model, with the aim of analysing the influence of quantum corre- 
lations on the spin-squeezing and entanglement effects. Section 5 summarizes the 
main results obtained in this paper and discuss some possible perspectives for fu- 
ture research. We have added three appendixes related to the calculational details of 
certain topics and expressions used in the previous sections: Appendix A shows how 
the unitary transformations associated with angular-momentum generators affect 
the expressions for variances, covariances, and RS uncertainty relation; Appendix 
B discuss the Kitagawa-Ueda model and the inherent spin-squeezing and entan- 
glement effects; while Appendix C exhibits an important set of specific numerical 
computations which leads us to establish a validity domain for the entanglement 
criteria studied in this work. 

2. Definitions and background for spin coherent states 

In this section, wc will briefly survey some mathematically appealing features inhe- 
rent to the su(2) Lie algebra with emphasis on a specific group of unitary transfor- 
mations involving the generators of this algebra which leads us to properly define 
the spin coherent states. For this task, let us initially consider the standard angular- 
momentum generators {3 x ,3 y ,3 z } which act on a finite-dimensional Hilbcrt space 
T-L N . The familiar commutation relation [J a , J&] = ie a bc.3 c for a, b, c = x, y, z (letting 
h = 1 for convenience), where e a b c denotes the Levi-Civita symbol associated with 
the three orthogonal directions, permits us to define the raising and lowering ope- 
rators through the particular decomposition J± := 3 X ± iJj, with [J + , J_] = 2J 2 
and [3 Z , J±] = ±J±. Moreover, J 2 := J 2 + J 2 + J 2 = J 2 + \ {J + , J_} characte- 
rizes the total spin operator and constitutes an important element of this algebra 
since J 2 (also known as a Casimir operator) commutes with all the generators 3 a . In 
particular, these results allow to construct, for example, a set of simultaneous eigen- 
states of J 2 and J z , that is, {\j,m) : \m\ < j , j — 0, \ : 1, |,.. .}, with well-known 
mathematical properties!^ 

(i) It is interesting to note that each particular member of this representation (here 
characterized by a specific eigenvalue j) is analogous to the Fock statd^of the 
electromagnetic field, since it can be created by the repeated application of the 
raising operator (remembering that 3 2 J. +1 = 0) on the 'vaccum state', 

_I 3 3 + m r\ 

\j, m) := %j 2 J+m j— — 77 1 j, - j) for %, s = . 

J J [j + my. s\(r — s)\ 
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Thus, the set m}} corresponds to a discrete, orthonormal and complete 
basis for the (2j + l)-dimensional vector space C 2j+1 of angular-momentum 
statesp^ whose completeness and orthonormality relations are expressed as 
follow: 

3 

^2 \j,m){j,m\ = l and (f ,m'\j,m) = S fJ 5 m ,. m . 

m=-j 

Consequently, the expansion of any quantum state belonging to this finite- 
dimensional space can now be prompty obtained, 

3 

I*) = c i )TO |j,m) 
m=—j 

where Cj. m = (j, m\^>) denotes the coefficients of such a discrete expansion, 
(ii) The nondiagonal matrix elements (J 1 , m'\ J± \j, m) associated with the moments 
J± show a dependence on the binomial coefficients ^ r . s , that is 

(j',m'\Ji\j,m) = f^L2±UlL\ 2 6j ,,S m/ , m±k (0 < k < 2j), 

while (J', m'\3^\j, m) = m k 8ji j5 m >, m demonstrates an explicit connection 
with the discrete component m. These results are extremely useful when ap- 
plied to the evaluation of nondiagonal matrix elements related to well-known 
families o f unit ary transformations involving the linear co mbination of spin 
operators J32E3] an( j a j so j n s t uc {y f S p m coherent states E 3 | 35 | 36 | 



2.1. Unitary transformations 

Let T(f2±,fi z ) :— exp (f2 + J + + fi z J z + S7_ J_ ) denote a kind of general abstract 
operator expressed in terms of generators of the su(2) Lie algebra and arbitrary 
c-number parameters fl± and Q z . For f2 + = £, f2_ = — £* and Q z = iu>, with £ € C 
and to £ K, such an abstract operator represents a generator of unitary transforma- 
tions whose respective generalized 'norma l'- and 'antinormal'-order decomposition 
formulae display the following expressions E31H1 

T(£, u) = exp (A+J+) exp [In (A.) J z ] exp (A_ J_) 

= exp (-A* J_) exp [- In (A*) J z ] exp (— A* J+) (1) 

with 

A (£/</>) sin(0) A (r/^)sin(0) 



cos(4>) — i (w/20) sin(</>) ' cos(0) — i (uj/2(p) sm(<p) ' 



,2 



A z =[cos(<p)-i(oj/2<j>)sm(<j>)] and </> = |£| 2 + (w/2) 

So, the action of T(£, cj) on the generators { J a } — here defined through the relation 
J a := Tt(£, w)J a T(£, ui) for a = x,y,z — produces effectively a new set of angular- 
momentum operators {J a } expressed in terms of the old ones that, by their turn, 
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are multiplied by determined coefficients which depend on the parameters £ and u>, 
namely, 

sin 2 (0)" 



cos(20) + 2Im 2 (£)^M 



•fx + 

sin(20) 



20 



-„f^ + 2 Re(£)Im(0^ 

„ Mf) i^ + ImK) 5i£M 



.... sin 2 (</>) „ sin(2(i) 
W lm(0— ^ +Re(0— ^ 



J,, 

J EC 



cos(2^) + 2Re 2 (0^|^ 



sin^) ^ sin(20) 
wRe(0— To MO — - — 



cos(20) 



w 2 sin (0) 



(2) 



(3) 



(4) 



Note that J 2 remains invariant under the unitary transformation ((T|), which implies 

in the identity J = J 2 ; besides, the geometric counterpart of this particular result 
is directly associated with arbitrary rotations on the surface of a sphere of radius 
+ Indeed, let T(0,ip) denote a particular case of T(£, w) when £ = | exp(— \ip) 
and u = 0, which implies in specific decomposition formulae characterized by A± 
■ tan(0/2) exp(=Fi(/3) and A 2 = sec 2 ((9/2). The connection between £ and A + reflects 
the stereographic projection of a two-dimensional sphere S 2 on the complex plane 
C with one-point compactification (in this situation, the infinite point corresponds 
to the north pole of S 2 = C U {oo}). 

As a consequence, the action of T(0, (f) on the vacuum state \j, —j) will define 
our next object of study: the spin coherent states — also known in current literature 
as SU(2) coherent states or spin coherent states!^ In fact, we will establish some 
few important mathematical results connected with the spin coherent states which 
constitute the first basic tools for discussing the subtle link between spin squeezing 
and entanglement. 



2.2. Spin coherent states 

In general, the spin coherent states can be defined by means of two different mathe- 
matical procedures: the first ond^ basically follows the Schwinger's prescription of 
angular momentum (note that LMG modeF^ represents a typical exam ple of this 
prescription) in order to produce such states, while the second on pursues an 
algebraic framework analogous to that adopted for the field coherent states which 
allows us to describe physical systems consisting of N two-level atomt0 (for instance, 

a See Refs. I42H47I for a detailed discussion on the atomic coherent-state representation and its 
practical applications in multitime-correlation functions and phase-space quasidistributions. 



I 



I 



November 27, 2012 2:12 WSPC/INSTRUCTION 



Spin squeezing and entanglement via finite- dimensional discrete phase space description 7 

see the Dicke modeF^j). Here, we follow closely the last way that consists in defining 
the spin coherent states through the mathematical relation \9, tp) := T(#, p)\j, —j), 
with the vacuum state \j, —j) written in terms of the previously mentioned basis 
m )}-j<m<j- Thus, after some algebra, it is easy to show that 

\0,<p) = [^2 3 , k sm 2k (9/2)co S 2 ^~ k \e/2)\ 2 exp(-ifcp)|j, k - j) (5) 

fe=0 

presents a bijective mapping with the TV-photon generalized binomial states related 
to the electromagnetic field, being this correspondence properly explored in Ref. 1491 
Next, we focus our attention on an important set of mathematical properties 
inherent to the spin coherent states ([5j that leads us to discuss the squeezing effects 
in determined physical systems under different circumstances. Note that a complete 
list of such properties can be promptly found in Refs. 1411 and [501 

Non-orthogonality. The inner product of two distinct spin coherent states can be 
evaluated through the completeness relation associated with {\j, m)}, yielding 
the well-known result 

(9',tp'\9, <p) = {cos(0') cos(0) + sin(0') sin(0) exp[i<y - p)]f' J . 

For 9' = 9 and ip 1 = ip, we achieve (0, tp\9, tp) = 1 (normalization condition). 
Now, excepting for the antipodal points, the spin c ohe rent states are in general 
not orthogonal. Hence, the Majorana-Bloch sphere^^ of unit radius represents 
the ideal geometric element to describe such states, since its respective north 
(9 = 7r) and south (9 = 0) poles correspond to the highest/lowest states \ j, ±j). 
For completeness sake, let us briefly mention that \(9' ,ip'\9,tp)\ 2 = cos 43 (9/2) 
(overlap probability) with cos(9) = cos(0') cos(6') + sm(0') sin(6>) cos(<^' — p) is 
limited to the closed interval [0, 1], which reinforces the overcomplete character 
of the states under consideration (they do not form an orthonormal set). 
Completeness relation. Let P(9,ip) = \9,<p)(9,<p\ denote the diagonal projector 
operator related to the spin coherent states which satisfies the property^ 

/ dfi(0, <p)P(6, <p) = l with dfi(0, p) = ?L±1 sin(6»)d6»d^. 
J ' 4?r 

In particular, such a completeness relation asserts that any quantum state |^) 
belonging to C 2:,+1 can be properly expanded in this overcomplete basis, namely 

l*> = J d£i{e,<p)*(9,tp)\e,<p). 

In this case, ^(9, <p) represents a polynomial function of degree 2j, whose ana- 
lytical expression is given by 

2j r n 1 
*(0, ip) = sin 2fc (0/2) cos 2 ^- k \9/2)\ 2 exp(ifc^) (j, k - 

fe=0 
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It is worth mentioning that operators acting on this Hilbert space also admit 
expansions in both the nondiagonal and diagonal forms, the diagonal represen- 
tation being of particular interest to deal with statistical operators for atoms.SH 
Minimum uncertainty states. Note that spin coherent states constitute an im- 
portant class of minimum uncertainty states. To dem onstrate this assertion, let 
us initially consider the RS uncertainty principld^ for 3 X and 3 y , i.e., 

- (r JxJ J 2 > \ j,]>| 2 = i|(j z )| 2 (6) 

where '^j x .j B = (\ {J x , Jy}) ~ (Jx)(Jy) represents the covariance function, and 
r ¥j a = (J 2 ,) — (J a ) 2 stands for usual variance when a = x,y. So, if one considers 
the previous mean values evaluated for the spin coherent states, it is immediate 
to verify the equality sign in this equation since both the expressions are equal 
to jj 2 cos 2 (8). In fact, this particular result can be extended in order to include 
any non-parametrized complex variable — see Appendix A for further results. 



3. Mappings via finite-dimensional discrete phase spaces 

In this section, we will introduce certain basic mathematical tools constituents of 
the pioneering formalism initially developed in Refs . IT5l and 1 1 71 for physical systems 
with finite-dimensional space of states W K . It is worth stressing that such tools will 
represent our guidelines for the fundamentals of the formal description of finite- 
dimensional phase spaces by means of discrete variables. For this specific task, let 
us first establish the mod(N)-invariant operator basifl^ 



v) := -j= ex P 



2vri 



N (w + M 



S(V,0 (7) 



expressed in t erm s of a discrete double Fourier transform of the symmetrized unitary 
operator basiJ^ 

Sfofl ^cxpf^W, 



where the labels rj and £ are associated with the dual momentum and coordinatelike 
variables of an N 2 -dimensional discrete pha se s pace here endorsed by an underlying 
presymplectic structure of geometric origin!^ Note that these discrete labels obey 
the arithmetic modulo N and assume integer values in the symmetric interval [—£, £] 
for £ = ^ fixed0 A comprehensive and useful compilation of results and proper- 
ties of the unitary operators U and V can be found in Ref. 1191 since the primary 
focus of our attention is the essential features exhibited by {G(n, i / )} J u ) i/=-^,...,^- 

b Henceforth, for convenience, we assume N odd throughout this paper. However, it is important 
to stress that even dimensionalities can also be dealt with simply by working on non-symmetrized 
intervals. 
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The set of N 2 operators {G(/i, v)} allows, for instance, to decompose any linear 
operator O acting on % K by means of the expansion 

l i 

°=N £ efa^Gfav), (8) 

where the coefficients 0(p,v) = Tr [GT(/i, i/)0] are evaluated through trace ope- 
ration and correspond, in this context, to a one-to-one mapping between operators 
and functions embedded in a finite phase-space characterized by the discrete vari- 
ables /i and i/PZl So, if one considers the density operator p in such a case, we verify 
that 

1 e 

P=K E ^v)G(^v) (9) 

admits a plausible expansion with intrinsic mathematical properties since the coef- 
ficients result in the discrete Wigncr function W p {p, v) := Tr [G^(p, v)p] . The first 
practical consequence of this decomposition yields the mean value 

1 1 

(O) := Tr[Op] = - £ 0{p,v)W p {p,v), (10) 

whose formal expression depends explicitly on the product of both mapped forms 
related to O and p. Consequently, the moments (J„) for a = x,y,z and < k < 2 j 
can now be promptly obtained from this formalism for any finite quantum state 
belonging to %2j+i - Next, we apply such a mapping technique for the spin coherent 
states in order to establish an important link between angular-momentum operators 
and (2j + l) 2 -dimensional discrete phase spaces. 



3.1. Applications 

According to a prescription adopted in Ref. ll8[ let {\j, m)}_j< TO <j denote the eigen- 
states of U with eigenvalues {oj m } for oj = exp ^^frr) nxe d- This assumption leads 
us to establish the general properties U' 7 \j, m) = u> mv \j, m) and \j, m) — \j, rn—^) 
(together with the relations U^' +1 = 1, V^' +1 = 1, and V?U r ' = w^UW e ), which 
represent a satisfactory mathematical connection for the present purpose. Next, we 
determine some useful results related to the angular-momentum operators and spin 
coherent states, remembering that G'(p, v) = G(p, v) for — £ < p, v < I. 
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As a first and pertinent application, let us consider the mapped expressions 
j 

m=—j 
j—k 

m=— j 

(3 h _){n,v)= J2 



(j + TO + fc)! f^2j,j+m+k 



(J + m)\ V ffzjj+m 

(j — m + k)\ ( J+m-k 



(j - to) 



m——j J r k 

which depend on the nondiagonal matrix elements 

27ri 



(j,m\G(ii,v)\j,m + k) 

^ ^^^^^^ 

(j,m|G(/i,i/)|j',m - fc), 
s * ' 

S^,m|f < m-*(/i.»') 



2j 



+ 1 ^ 



exp 



2j + l 



— to') + (m' — to) ( ^ + — 



associated with the discrete mapping kernel G(/i, i/) for m — m±k and < k < 2j. 
Since the diagonal matrix element m u m (/*, v) — $m^^ gives the Kronecker delta 
function (in such a case, the superscript [2j + l] denotes that this function is different 
from zero when its labels are congruent modulo 2j + 1), it seems reasonable to 
obtain a simple expression in the first example, that is (J*) (/•*, v) — /-t ; however, 
this fact is not verified for (J±) (/U, v~). From the kinematical and dynamical point 
of view, these mapped expressions represent an advantageous set of mathematical 
tools which allows us to co mpre hend certain intriguing problems related to discrete 
symmetries of spin systems 

Basically, the second application consists in evaluating the discrete Wigner func- 
tion ## i¥ ,(/x, v ) present in the expansion of the projector P(0, (p) via Eq. (|9]), namely, 



In general, this particular quasidistribution can be written in terms of its respective 
discrete Weyl function Wo,tp(r],£) ■= T¥[S(r7, £)P(#, <p)] as follows: 



V2j 



1 j 

= y 



exp 



2?ri 
2j + l 



(ii) 



So, after some lengthy calculations, the analytical expression for We^ir}, £) assumes 
the exact form 



Av,0 = 



i 



V2J 



= J2 ex p 



2vri 



?7 TO 

2j + 1 ' V 2 



(6,<p\j,m- £,)(j,fn\9,(f), 



where (j, m\6,ip) can be promptly obtained from Eq. (|5|). It is interesting to stress 
that both the north (9 = tt) and south (9 = 0) poles of the Majorana-Bloch sphere 
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correspond to the limit cases W Wt¥> (u, v) = 5 3 - ' and Wq. v {}i,v) — , which 

reflect directly the highest/lowest states \j, ±j) for any ip £ [0,2ir). Furthermore, 
numerical calculations related to Eq. (|11[) suggest that some relevant contributions 
are connected with discrete values of the angles 6 and ip — in particular, integer 
multiples of 2 j+i • This evidence supports the ideas of Buniy and coworkers^ about 
a 'possible discretization' of the Majorana-Bloch sphere, as well as reveals some 
important features — and not yet properly explored — on the discrete nature of 
quantum states in finite-dimensional Hilbert spaces]^ 



3.2. Time evolution 

Here, we establish a mathematical recipe that permits to investigate the dynamics 
of a particular quantum system characterized by a finite space of (discrete) states. 
For this task, let p{t) describe the state of this physical system whose interaction 
with any dissipative environment is, in principle, automatically discarded. Besides, 
let us consider (for convenience) only time-independent Hamiltonians; consequently, 
the time-evolution of that density operator will be governed, in such a case, by the 
well-known von Neumann-Liouville equation ihdtpit) = [H, pit)]. So, those initial 
premises represent the constituting blocks for the aforementioned recipe. 

As a first proposal in our prescription, we determine a mapped expression for the 
von Neumann-Liouville equation which describes the time evolution of the discrete 
Wigner function W p (u, v;€). In this sense, the mapping technique sketched in this 
section yields the differential equation 

t 

m t w p {^v-t) = ^h(^,mVTOmV;0, (12) 

li' ,v'=-t 

where Jz? h (m, v, A*'i v ') represents the mapped form of the Liouville operator written 
in terms of H (//', v") = Tr [G(/x", v")U], that is, 

Sf a (jM, v, fx', v') = -^ J2 sin [^("Z 3 ' " a 'P)\ ex P { ^ W» ~ /"') + P{v ~ v 1 )} \ 

x exp |^ [a'(p - (i") + - v")}} H{jx", v"). 

Note that A denotes the set {a, (3, a 1 , /?', //', v"} £ [—£, £} in this expression. In the 
following, let us mention some few words about the formal solution of Eq. (fT2l) : it 
can be expressed analytically in terms of the series^ 1 

i 

W p (n,is;t)= J2 ^(^^t\K,T;t )W p {K,T;t ), (13) 



c In the quantum-gravity scope(531it is worth mentioning that certain theoretical approaches to the 
generalized uncertainty prin cipl e (GUP) also suggest the breakdown of the spacetime continuum 
picture near to Planck scalel^l From an experimental point of view, it seems reasonable to argue 
that 'discreteness is actually less speculative than absolute continuity'. 
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whose N 2 -dimensional discrete phase-space propagator admits the expansion 
3>(M, v- t\K, r; t ) = 8^8™, + ^(t - *o)JSf H (/i, k, r) 

+ ^(^-<o) 2 51 -^h(m,^<t')^h«t' ! k,t) + --- 

K',T'=—t 

which allows to evaluate directly the time evolution of v\ t) by using the series 

related to the iterated application of the mapped Liouville operator. The advantages 
and/or disadvantages from this particular formal solution were adequately discussed 
in Ref. 1201 and subsequently applied with great success in the discrete Husimi- 
function context for the LMG model!^ 

The alternative proposal establishes a differential equation for the discrete Weyl 
function analogous to Eq. (TT2"j) . Then, after some calculations related to the product 
of three symmetrized unitary operator bases and its respective trace operational^ 
it is immediate to show that such an equation can be written as 

e 

md t w p ( v ,tt)= ]T £ H (^,r/,£')W, (14) 

n',£'=-t. 

where W p (r), := Tr [S(rj, £,)p(t)] corresponds to our object of study and 



2i 



N 



N 



Tr[S(>j-»j',£-£')H] 

yields the mapped forms of the Liouville and Hamiltonian operators. As expected, 
Eq. (|14p describes the time evolution of the discrete Weyl function whose formal 
solution is given by 

e 

ir p (v,t;t)= £ l P(v,^t\v'^';toW P (r 1 '^';t ), (15) 

with V(r), f; t\rf , to) being the discrete dual phase-space propagator which admits 
a time expansion similar to &(fJ>, v\ t\n, r; to), that is 

V(r,, £; t\ V ', to) = 8%]/^ + ^(t - i Q )£ H (r?, £, rf , O 

■2 t 

+ 2w {t ~ to)2 E £*(v,z,v",e')£*(v",e,v',t;')+--- ■ 

v",i"=-t 



Since the discrete Wigner and Weyl functions are connected by means of a double 
Fourier transform, the low operational costs involved in this mathematical recipe 
are advantageous — from a computational point of view — if compared with the 
previous one. Notwithstanding the apparent advantage, it is worth remembering 
that W p {rj,^;t) is a complex function and hence it represents an intermediate step 
for the main goal, namely, the time evolution of the discrete Wigner function — and 
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consequently, the time evolution of certain mean values via Eq. (|10p . The schematic 
diagram shown below 



FT 



represents a summary of the previous proposals for describing the dynamics of a 
finite quantum system represented in a N 2 -dimensional discrete phase spacePSl 



4. Spin squeezing and entanglement 

This section will illustrate how the mapping techniques leading to finite-dimensional 
discrete phase spaces can be effectively used in the study of important quantum 
effects, such as spin squeezing and entanglement. Basically, we will invoke the LMG 
model already discussed in current literature (e.g., see Ref. |40| which exhibits not 

IOQI 

only strong mathematical and physical appeals}^ 1 but also a fundamental physical 
property: quantum correlations due to different types of interactions and mediated 
by a transverse magnetic field. Since quantum correlations are responsible for the 
aforementioned effects, it is natural to investigate how the interactions involved in 
this model affect the standard quantum noise related to the spin coherent statesEI 



4.1. The modified Lipkin-Meshkov-Glick model 

Here, we adopt the prescription established by Vidal and coworkers^ 1 for the modi- 
fied LMG model through the Hamiltonian operator (written in the spin language) 



H = -2hJ z - 2p H 



(N/2)l -p-fJi+Ji), (16) 



which describes a set of N spins half mutually interacting in the xy plane subjected 
to a transverse magnetic field h. The coefficients p± = (I1L7) have an important 
role in this description: they allow to (i) investigate the different anti-ferromagnetic 
(A < 0) and ferromagnetic (A > 0) cases inherent to the model for any anisotropy 
parameter I7I < 1 (7 = 1 refers to the isotropic case), and also (ii) ensure that the 
free energy per spin is finite in the thcrmodynamical limit. In particular, such a 



d I n A ppendix B, we investigate a soluble spin model with lowest-order nonlinear interaction in 
J z (SI whose algebraic features permit to attain a particular set of analytical res ults very useful 
in the study of spin squeezing and entanglement effects via quantum correlations PJJ 
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model presents a second-order quantum phase transition at A = for A > fixed, 
whose symmetric (A < \h\) and broken (A > \h\) phases are well-defined within the 
mean-field approach; furthermore, its ground-state entanglement properties exhibit 
a rich str ucture which reflects the internal symmetries of the Hamiltonian operator 
H.EHEH Indeed, such an operator preserves the magnitude of the total spin operator 
since [H.J 2 ] = for all 7, and does not couple states having a different number of 
spins pointing in the field direction, namely, 

N -, 

H, J^J erj 2 =0 (spin- flip symmetry). 



Consequently, it is immediate to verify that Eq. ([16]) can be diagonalized within 
each (2j + l)-dimcnsional multiplet labelled by the eigenvalues of J 2 and J^, this 
fact being responsible for the soluble character of the associated spin modeljj 

Next, let us mention some few words about two discrete conserved quantities 
inherent to the spin model under investigation which reflect certain additional sym- 
metry properties for 7 = 1 fixed. In this case, the simplest operator commuting with 
H, therefore giving a constant of motion, is the parity operator II — here defined 
as II := R z (7r) = exp(i7rJ 2 ). This result tells us that the Hamiltonian matrix, in 
the J 2 representation, breaks into two disjoint blocks involving only even and odd 
eigenvalues of J z , respectively. The second interesting quantity comes from the anti- 
commutation relation {H, R} = for R(7r, 0, 7r/2) = H x (tt)IL z (tt/2): it corresponds 
to a particular rotation of the angular-momentum quantization frame by the Euler 
angles (7r,0,7r/2), transforming, in this way, H — s- H for the isotropic case. Thus, 
if \Ej) is an energy eigenstate with eigenvalue Ej, then R\Ej) is also an eigenstate 
of H with eigenvalue — Ej. This specific symmet ry property gives rise to an energy 
spectrum that is symmetric about zero PAUSES] 

After this condensed review, we establish below a sequence of steps that allows 
us to calculate the time evolution of several quantities necessary for investigating 
both the squeezing and entanglement effects associated with the aforementioned 
spin model in terms of the parameters h and p±. The first one consists in adopting 
the theoretical framework described in the previous section for the time-dependent 
discrete Wigner function defined upon a (2j + l) 2 -dimensional discrete 

phase space labeled by the angular-momentum and angle pair (/x, v) € [—j, j] — in 
particular, we adopt the theoretical prescription established in Ref. [20] for the angle 
variable, that is Q v = v. Since this approach depends on the Wigner function 
evaluated at time to = 0, the second step consists in fixing the spin coherent states 



c It is important to stress that the prescritpion here adopted does not coincide with the origi- 
nal Li pkin -Meshkov-Glick model^H] which was initially introduced over 40 years ago in nuclear 
physics®S for treating certain fermionic systems. In fact, this new modified version brings to scene 
the modern statistical mechanics point of view, where the collective properties of spin systems can 
be worked out with great success. ^J>5] The extended original version of the LMG model, which 
explores the parity symmetry via finite-dimensional discrete phase spaces, can be found in Ref. 1181 
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as initial state — see Eq. (fTTj) . The last step refers to the numerical calculation of 
the moments 

(JSWW = E (17) 

n,v=-j 

and covariance functions in the Schrodinger picture fixing the initial state in 9 = \ 
and ip = (in this way, the uncertainties are redistributed between the orthogonal 
components in the yz-plane), and also considering (for convenience, not necessity) 
only the ferromagnetic case A = 1. Henceforth, Eq. (fl6|) will assume the simplified 
form H' = —hJ z — jj (j^ + 7J?), where the constant term |(1+7)I was suppressed 
at this initial stage since it represents only a phase for the time evolution operator 
with null contribution (as expected) within our computational approach. 

4.2. The match between squeezing and entanglement effects 

In this subsection, we adopt the theoretical framework established for the Kitagawa- 
Ueda model (see Appendix B) concerning the match between squeezing and entan- 
glement effects. In particular, let us initially consider the results obtained in Table 
[I] for the RS uncertainty principle related to the angular- momentum generators and 
its connections with the 5-inequalities. Since the covariance function i / j a j b (t) has a 
central role in this algebraic approach and requires only numerical computations of 
(^a)s,ip(t) an d ({Ja> ^b})e, ip{t)i it is convenient to clarify certain subtle steps that 
are inherent to the exact analytical calculation of the last quantity. According to 
the mapping technique discussed in Section 3, the expression for 

1 3 

({Ja,J6}>^(*) = ^t-j J! {{*a,3b})(wWe, v {u,v-t) (18) 

u,v=-j 

basically depends on the mapped form of the anticommutation relation, that is 
({J Q , J b }) (/,, v) = 1 £ r A ( M , u\a', 1/, n", v") (J B ) (/,', v 1 ) (J b ) (/i", v") 



where Q stands for {//, 1/ , fi", z/'} G [—j, j] . It is worth stressing thatP^ 
Ta(p, v\h', v', n", v") = — ^ ^ 2 Y2 ex P j 2~^Y - MO + CO ~ 



— l>i^[ri"(»-»")+Z"(v-v")} 



2j + l 

shows explicitly the embryonic structure of the continuous cosine function presents 
in the well-known Weyl-Wigner-Moyal phase space approach!^ — here, Q' denotes 
the set {?/,£', ?/',£"} € [— Thus, the iS-inequalities can be properly estimated 
for the modified LMG model, which lead us to investigate the squeezing effects. For 
completeness sake, let us briefly mention that ([3 a ,^b]}e,tp(t) yields an expression 
analogous to Eq. (|18[) , but with two minor modifications: the function cos(z) in Ta 
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Table 1. The 5-inequalities exhibit unique mathematical virtues since they yield a direct con- 
nection with the RS uncertainty principle related to the non-commuting pair {J a , Jb}a,b=x,y,z 
of angular-momentum operators, where the covariance functions have an important role. Indeed, 
such inequalities lead us to investigate the squeezing effects through a well-established criterion 
in literature: for S^' < 1 (squeezing condition) and S^ c ' > 1, the inequality S^S^ > 1 is 

( c ) —(c) 

always preserved; moreover, the saturation S a = 1 describes minimum uncertainty states. 

In this table, we show all the possible links among RS uncertainty principles and cS-inequalities, 
with "fj a j b restricted to the closed interval [— ^/"l / 3 a "f / j b , ^/"^j a "l / j b ] ; moreover, the superscript 

(c) of the product S a ' denotes the angular-momentum component resulting from the com- 
mutation relation between J a and 3 a . 



RS uncertainty principles 

{Jz, Jy, J 2 } £= %N 




■^.-denominators 

K abc + o Vp e w N 




^-inequalities 

s^sf ■.= r 3 .-rj t /ii' abe 


Vj„-(*W 2 > *l(J«>l a 

*.%,-(%.J.) 3 > jl<j»>l a 

Vj* -(%,j.) a > 2K J ->I 2 


t^xyz 
t^xzy 
t^-yzx 


:= [(*W 2 + 3 K J *>I 2 ] 
== [(Ti.J.) a + ||(J,>| a ; 

■= [W B JJ 2 + |K J -)I 2 ] 


1 

2 

1 

2 

1 
2 


s^s« > 1 
Sx y) si y> > 1 

Sy^S^ > 1 



should be replaced by i sin(z) (keeping constant the argument of the trigonometric 
functions), which implies in the change Ta — > Tc- 

In order to illustrate the theoretical framework developed in this work, let us now 
consider the 5-inequality SySi x ^ > 1 and the ^-inequalities {£ y ,£ z } described in 
Appendix B by means of Eq. (|B.3[) . In particular, this procedure allows not only to 
establish a direct link between squeezing and entanglement effects for the modified 
LMG model, but also to investigate how these effects are affected by the transverse 
magnetic field h and anisotropy parameter 7 with N = 20 fixed. For instance, 
Fig. EJa,c) represents the plots of S y x \t) (dot-dashed line) and Si (t) (solid line) 
versus t G [0,50] for different values of (^,7): (a) (-0.1,0.2) and (c) (-0.13,0.1). 
It is interesting to observe at a first glance how the squeezing effect is sensitive 
to small variations in the parameters \h\ and 7. Further numerical computations 
corroborate this sensitivity and allow us to produce the following proper description 
about such an evidence as time goes on: for 7 > \h\, both the parameters Sy (t) and 
S {x) {t) exhibit squeezing effect and oscillatory behaviour due to the prevalence of 
quantum correlation effects (here introduced by the anisotropy parameter 7) on the 
transverse magnetic field h; however, if one considers 7 < \h\, the squeezing effect 
is strongly reduced (in fact, it always survives for small values of time) in both 
the parameters, which reveals the relative strong influence of h on the quantum 
correlation effects. With this in mind, let us now consider Fig. [2jb,d) where the 
plots of £ y (t) (dot-dashed line) and £ z (t) (solid line) are exhibited for the same 
values of (h, 7) used, respectively, in[^a,c). The 'almost perfect match' between 
squeezing and entanglement effects demonstrate the previous assertions and indeed 
reinforce the fundamental sequence correlation 1— > entanglement t— > squeezing of 
quantum effects. 

At this moment, we pay special attention to the results obtained by Toth and 
co-workers^ where a full set of generalized spin-squeezing inequalities for entan- 
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Fig. 2. The 'almost perfect match' between squeezing and entanglement effects is clearly visible 
in this selection of pictures, where the parameters (a,c) Sy X ^ (t) (dot-dashed line) and <s£ (t) (solid 
line) as well as (b,d) £ y (t) (dot-dashed line) and E z (t) (solid line) are plotted as functions of time 
and placed side by side for immediate comparision. Furthermore, this set of plots corresponds 
to two different values of (h, y), namely, (a,b) (—0.1,0.2) and (c,d) (—0.13,0.1), which leads us 
to investigate how the quantum correlations are affected by the transverse magnetic field h and 
anisotropy parameter 7 via (2j + l) 2 -dimensional discrete phase space inherent to the modified 
LMG model. 



glement detection was derived and studied in detail. For such task, let us formally 
establish an entanglement criterion for separable states through a set of inequali- 
ties involving the mean values { (J a ), (J„) } a - x y z that summarizes those aforemen- 
tioned results. 



Entanglement criterion. Let p describe an iV-particle density operator, as well 
as {(Ja;), (Jj/), (J^), (J|), (Jy)> (Jz)} characterize the mean values of collective 
angular-momentum operators associated with the physical system of interest 
which, by hypothesis, are previously known. For separable states represented 
by the convex sum p — J2k PkPk ® Pk ® • • -® Pk N ~^ ® Pk^ > wnere Pk features 
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the probability distribution for a given k, one verifies that 

<J*> + <J*> + (J*) < \N(N + 2) (19) 

i (N + 2) (r 3m + % y + r Js ) > JiV(iV + 2) (20) 

|(JV + 2) [{3 2 a ) + (J 2 h ) -(N- l)Y 3c ] < ±N(N + 2) (21) 

(N — 1) (r Jo + r Jb ) - <J2) + JV > ^N(N + 2) (22) 

are always fulfilled, \N(N + 2) being considered a lower/upper bound in all 
cases, with {a, 6, c} labelling all the possible permutations of {x,y, z}. Hence, 
the violation of any inequality (|2"0"|) - (|2"21 leads, in this framework, to entangled 
states since Eq. (fT9|) is valid for all quantum states — indeed, such inequali- 
ties define a polytope in the three-dimensional ((J^), (Jy), (Jz))-space, where 
the separable states lie inside this geometric representation. It is important to 
stress that physical systems with particle- number fluctuations (i.e., the particle 
number N is not constant, such as in BEC experiments) are discarded in this 
criterionp3 

Remark 1. Let us now apply these inequalities to the modified LMG model with 
the main aim of corroborating the previous results and presenting alternative forms 
for detection of spin-squeezing and entanglement effects. Numerical computations 
show that Eqs. ([21]) and (|22|) are not satisfied for determined values of time, which 
leads us to choose one of them whose particular time evolution combines with those 
exhibited by S and £ -parameters — see Fig.[2][a)-(d). So, let us define the parameter 

. (n - i)r Jo 

E a ■= /T2X ~ /T2X _ jv > 1 <"' r = .'/■ '- ) (2-->) 



(J 



based on Eq. (|2Tj) . where the violation E a < 1 directly implies in the quantum 
effects under investigation. Note that, compared with Eq. (|B.3|) . E a has subtle 
differences: it requires the sum of mean values related to the angular-momentum- 
squared operators (see denominator of both expressions), while the previous one 
involves only squared mean values of the angular-momentum operators. Figure [3] 
shows the plots of E y (t) (dot-dashed line) and E z (t) (solid line) versus t G [0, 50] for 
the same set of values associated with the transverse magnetic field and anisotropy 
parameters used in the previous figure, i.e., (a) (—0.1,0.2) and (b) (—0.13,0.1). It 
is important to stress that the functional characteristics of the involved parameters 
in the study of spin-squeezing and entanglement effects, as well as the associated 
entanglement dynamics to the modified LMG model, represent two essential factors 
that could explain (since the time evolution also depends on the initial state chosen 
for the physical system), in part, the similarities of behaviour related to the S, £, 
and .©-parameters viewed in Figs. and [3J In summary, the entanglement criterion 
conceived by Toth et aZ.EGIl reinforces the 'almost perfect match' between spin- 
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(a) (b) 




Time Time 



Fig. 3. Plots of E y (t) (dot-dashed line) and E z (t) (solid line) as a function of time for the same 
values of (/i,7) used in the previous figure, that is, (a) (—0.1,0.2) and (b) (—0.13,0.1). Here, two 
factors immediately emerge from our considerations on these pictures: (i) their similarities with the 
previous cases (emphasizing, mainly, the strong influence of h on the entanglement effects), and 
consequently, (ii) the important role of quantum correlations in such measurements. It is worth 
mentioning that the particle number N used in the numerical calculations remains the same (i.e., 
N = 20), which preserves important basic quantum effects (contrary to the thermodynamic limit). 



squeezing and entanglement effects, which leads us to investigate, in this particular 
case, the quantum correlation rules via discrete Wigner and Husimi functions@ 

Henceforth, let us mention some few words on the entanglement dynamics of the 
Hamiltonian operator H' (Schrodinger picture): it basically involves the action of 
collective spin operators responsible for introducing the quantum effects attributed 
to two-body correlations — here mediated by the anisotropy parameter 7 — upon a 
nonseparable initial state, that is, the collective spin-coherent state ([5|) with 9 = ? 
and ip = fixed. In this sense, since the discrete Wigner function reflects the action 
of the time-evolution operator U(i) := exp(— -^H't) upon the initial state |f,0) 
in a (2j + l) 2 -dimensional discrete phase space, it seems natural to investigate its 
behaviour for those specific values of time where effectively occur the spin-squeezing 
and entanglement effects. 



4.3. Finite- dimensional discrete phase spaces: a case study of 
quantum correlations via Wigner and Husimi functions 

Initially, let us consider the time evolution of v; t) upon a 21 2 -dimensional 

phase space labeled by the pair v) of dimcnsionlcss discrete angular-momentum 
and angle variables where the particular operator H' = — hJ z — (J 2 , + 7J 2 ,) will 
take place at this description. As usual, the parameters employed in the numerical 
computations are the same as those mentioned in figure [^b), i.e., N = 20 and 
(h, 7) = (—0.1, 0.2). Furthermore, we choose six representative time values of £ z {t) 
which reflect the first local minimum and maximum points exhibited in figure J^b). 

f For completeness sake, we present a case study for h = (absence of transverse magnetic field) 
in appendix C, whose numerical computations lead us to discuss the real necessity of demanding 
the respective validity domains of the entanglement criteria used in this work. 
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5 10 15 20 5 10 15 20 5 10 15 20 



Fig. 4. (Color online) Time evolution of Wz. v;t) for the modified LMG model with N = 20 
and (h, 7) = (—0.1, 0.2) fixed, where the labels fi and v characterize, respectively, the dimensionless 
angular-momentum and angle pair. In particular, these pictures show how the two-body interaction 
term present in the Hamiltonian operator H' affects the initial Wigner distribution q(h,v;0) 
(via contour plots) for determined values of time. We have adopted in our numerical computations 
the values (a) t = 0, (b) t = 2.15, (c) t = 4.75, (d) t = 7.10, (e) t = 9.05, and (f) t = 9.95, which 
illustrate a representative but not complete evolution of the discrete Wigner function. Moreover, 
these values correspond to the first local minimum and maximum points of £ z (t) — see solid 
curve in figure[2fb) — while (f) represents the situation £ z (9.95) ~ 1. Although this phase space 
appears to the eyes as continuous (due to the large number of points associated with N = 20), it 
is important to stress that such a particular phase space is (by construction) genuinely discrete. 

In this way, figure [4] shows the contour plots of W* t o([i, u; t) versus (/x, v) 6 [—10, 10] 
where, in particular, (a) t = (initial Wigner function), (b) t = 2.15 (first local 
minimum point), (c) t — 4.75 (first local maximum point), (d) t = 7.10 (second 
local minimum point), (e) t = 9.05 (second local maximum point), and (f) t = 9.95 
(corresponds to £ z (9.95) ss 1). In all these pictures, it is interesting to observe how 
quantum correlations associated with two-body interaction term modify the initial 
correlations present in the collective spin-coherent state |f,0); besides, note that 
both fi and v were conveniently shifted in the snapshots. Next, let us describe certain 
specific points inherent to the time evolution of the discrete Wigner function. 

• 4(a). The function Wjl^{^, 1 V] 0) exhibits a rotational symmetry with negative 
values located at the proximities of (—10, 0) and (10, 0) (see small orange zones); 
as well as different widths in the respective fi and ^-directions. This important 
fact suggests, if one considers both the directions, that marginal distributions 
related to the angular momentum and angle variables can be somehow relevant 
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and necessary in the analysis of quantum correlations. The approximate format 
of an ellipse — white zone in the middle part of the contour plot — represents a 
region where there exists raised probability peaks which feature, by their turn, 
the initial correlations associated with | f , 0) . 

• 4(b). This first case of time evolution corresponds to l W^fi{\i,v\t = 2.15) and 
exemplifies the first local minimum point of £ z (t) (see solid curve in figure 
[2th) when t = 2.15). The entanglement effects (associated with the two-body 
interaction, and synonymous of quantum correlation) introduce additional sym- 
metries that modify the statistical weights of the initial quantum state | § , 0) , 
originating, in this way, new interference patterns which should explain the 
changed widths of the discrete Wigner distribution function in t = 2.15, as well 
as the appearance of new zones where negative probabilities occur. Note that 

v; t) exhibits a motion towards the frontier of the discrete phase space 
labeled by fi = —10, the transverse magnetic field h being the agent responsible 
for such behaviour. Numerical computations indeed corroborate this assertion 
and allow us to show that, for h = 0, the Wigner function remains frozen in 
the discrete phase space, having its widths modified as time goes by. 

• 4(c). This particular case depicts ^^(/U, v\t — 4.75) and represents the first 
local maximum point of £ z {i). Here, W* t o(f/,, V] t) begins its motion towards the 
centre of such a discrete phase space, where the interference patterns associated 
with the components of XJ(t — 4.75)|f , 0) have specifically yielded the portrait 
verified in this situation. The small white 'island' located in the right hand side 
of this picture and quite near to the raised probability peaks (main white zone) 
depicts a small region with negative probabilities. 

• 4(d). The second minimum point of £ z (t) has as discrete phase-space repre- 
sentative the contour plot of W^.fi{^i,v;t — 7.10). These unique interference 
patterns (once H' describes for h ^ a nonperiodic system) with small 'is- 
lands' of positive and negative probabilities represent the accumulated effects 
attributed to the transverse magnetic field and anisotropy parameter. Moreover, 
let us briefly mention that XJ(t = 7.10) |5, 0) has a coincidence probability (also 
known as fidelity and/or time correlation) with the original state close to 89%, 
while BJb) results in 55%. 

• 4(e). Note that W^fi{fX, u;t = 9.05) corresponds to the second local maximum 
point of £ z (t), he-, we are describing a situation where the entanglement does 
not occur. For the sake of comparison, the coincidence probabilities associated 
with IDJc) and E£e) assume the respective percentages of 48% and 78%. This 
difference can be justified, in principle, by means of a simple inspection between 
both the pictures: while |4je) is located near the centre of discrete phase space, 
IDJc) stays in the proximity of /i = —10. However, if one considers the previous 
case [UJd) for subsequent comparison, different contributions should be taken 
into account in the analysis. 

• 4(f). This last case depicts W^^iji, v;t = 9.95) for £ z (9.95) « 1. It is interesting 
to observe how the nonperiodic dynamics related to the Hamiltonian operator 
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H' affects any reconstruction process of #^..o(p, v\ 0): some numerical estimates 
of $(t) ■= Tr[p(0)p(t)] for t — 9.95 results in 51% of time correlation. 

Now, let us show an essential formal result that explores the connection between 
discrete Husimi and Wigncr distribution functions}^ 

1 3 

^,i/;t) = __ I j2 £(/i,*vy) wx;*)- (24) 

Here, E(/j,, v\a ',y') defines a smoothing process characterized by the discrete phase- 
space function^ 

that closely resembles the role of a usual Gaussian function in the continuous phase 
space, since Jf(i], £) := £)/^#(0, 0) is here responsible for the sum of products 

of Jacobi ^-functions evaluated at integer arguments, namely, 

= ^{0a(«?M [0 3 («*I«O + expMM^M] 

+ exp(i7r0^4(ar?|ia) [#3(a£M + exp[i7r(?7 + 2j + l)]# 4 (a£|ia)]} 

where a = (4j + 2) _1 — in particular, see Appendix A of Ref.UHlfor definitions and 
technical details on Jacobi ^-functions. So, in order to illustrate how the smoothing 
process acts on the summand of Eq. (j2"l|) . let us consider #^l i0 (/^> v \ t) and also the 
same set of parameters used in the previous figure. In addition, it should be noticed 
that Mkfi(fi, v; t) is strictly positive as well as limited to the closed interval [0, 1] for 
any t > 0; consequently, all the negative values and non-regular patterns described 
in figurc[4]will be washed out by E(fi, v\fi', v'). Figure [5] exhibits the contour plots of 
Jifk t o(fj,, v\ t), where such effects are effectively checked along the pictures [5ja-f). In 
these cases, the width changes observed for the (/x, ^-directions are related to the 
two-body interaction term (here mediated by means of N and 7), while the motion 
towards frontier located at p = —10 is associated with the transverse magnetic field 
h. Since the widths associated with the angular-momentum and angle distributions 
are also affected by those effects, it seems convenient at this moment to employ the 
Wchrl-type entropy functionala 1 ^ for understanding how the correlations between 
the discrete variables fj, and v of a finite-dimensional phase space change for t > 0. 

To this end, let us now consider the mutual correlation functional I[H; t] defined 
as follows!^ I[H; t] := E[Q; t] + E[R; t] - E[H; t] > 0. In this particular definition, 

1 3 - 

E[H;t] = - — 1 ^M^W-Kit*,^)} (25) 
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Fig. 5. (Color online) Time evolution of o(m, v\ t) for the same set of parameters employed in 
the previous figure. In what concerns the discrete Husimi function, it must be stressed that such 
a function is particularly obtained from Ws. o(fj,,v;t) through a smoothing process characterized 
by a discrete phase-space function that closely resembles the role of the Gaussian function in the 
continuous phase space. This process essentially explains the disappearance of all irregular regions 
and negative values exhibited by the discrete Wigner distribution, which implies in a well-behaved 
function along its time evolution and limited to the closed interval [0, 1]. Similarly to the previous 
case, it is worth mentioning that such a particular phase space is genuinely discrete (by ab initio 
construction), even though it appears to the eyes as contionuous. 

corresponds to the time-dependent joint entropy functional here expressed in terms 
of the discrete Husimi function (l24l). while 



(26) 



1 3 

E[Q;i] = --_= ^(m;*)M^p(m;*)] 



and 



1 3 

E[R;i] = = J2 # P (r,t)hi[a fi (v;t)] 



(27) 



represent the marginal entropies — which are directly related, by their turn, to the 
respective marginal distribution functions (for technical details, see Ref. [2Tj) 



1 J 1 J- 

r^j M 'p^v\t) and 3$ p {v;t) = , = ^ Jt? p (fj,,v;t) 



V2j . .. 

y=-j i*=-3 

Note that E[H;t], E[Q;t] and E[R;t] constitute some basic mathematical elements 
for describing functional correlations between fi and v, and consequently, the width 
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Fig. 6. Plots of (a) E[H; t] and (b) I[H; t] versus t £ [0, 50] for two different values of transverse 
magnetic field h and anisotropy parameter 7 with N = 20 fixed, namely, the solid curves in 
both pictures correspond to (—0.1,0.2), while the dot-dashed curves are related to (—0.13,0.1). 
In each case, the orange region describes particular situations of JPm. o(fJ., v\ t) where the entropy 
functional are restricted to (a) E[H; t] < E[H;0] RJ 0.1994 and (b) I[H;t] < I[H;0] ss 0.7144. It is 
important to stress that E[H;t] and I[H;t] are sensitive to small variations of h. 



changes associated with the discrete Husimi function Jf p (fi, v\ t). Within several im- 
portant properties inherent to the entropy functionals, the Araki-Lieb inequality^! 
|E[Q;i] - E[R; t]| < E[H;i] < E[Q; t]+E[R; t] has a central role in such a description: 
for instance, it leads us to investigate how the dynamic correlations (introduced by 
means of the Hamiltonian operator H') affect the underlying correlations of the ini- 
tial state p(to) — here mapped onto J4? p (p, v\ to). So, if one computes the entropy 
functionals from J&%fl(p,i/',t), some interesting results can be promptly achieved. 

Figure[S]illustrates the time evolution of (a) E[H; t] and (b) I[H; t] as a function of 
t € [0, 50] for the same values of (h, 7) used in figure[2] Then, let us initially focus on 
(—0.1, 0.2) (see solid curves in both the pictures) since this case coincides with that 
used in the contour plots of J^^Qi, v\ t) for specific times. Those concentric ellipses 
observed in figure [5ja) can be considered as a phase-space signature of the initial 
state |§,0) at t = 0, here endorsed by E[H; 0] « 0.1994 and I[H; 0] « 0.7144, which 
means that such a state presents an initial correlation between \i and v featured by 
a spread angular-momentum distribution and 'localized' angle distribution. Indeed, 
the sum E[Q; 0] +E[R; 0] sa 0.5150 justifies such an assertion and (as expected) also 
corroborates the right-hand side of the Araki-Lieb inequality. Already at t = 2.15, 
one observes the displacement of the discrete Husimi function towards the frontier 
H = — 10 accompanied, in this situation, by a decrease of asymmetry between its 
intrinsic widths in the (/x, z/)-directions. Since E[H; t — 2.15] ~ 0.1958 implies in ef- 
fective information gain upon the discrete angular-momentum and angle collective 
variables, the extra correlations introduced through the two-body interaction term 
explain, in principle, the effects depicted in pictures (5Ja,b). However, the deformed 
widths viewed in(5][c-f) after the first rebound from that frontier and subsequent 
return towards the centre of the finite-dimensional discrete phase space are respon- 
sible for the inscreased values of both entropy functionals when 4.75 < t < 9.95. 
Let us now briefly discuss the case (—0.13,0.1) here characterized by a significant 
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reduction of the effects related to 7. Although the information gain happens more 

frequently in such a case, the spin squeezing and entanglement effects exhibited in 

pictures EJc,d) are associated with frequent motions of J4?kfi(fj,,L';t) towards the 

frontier fi = —10, followed by small variations of its widths. 

A last pertinent question then emerges from our considerations about discrete 

Husimi function and Wehrl-type entropy functionals: "How the smoothing process 

characterized by Eq. (f24|) affects the description of correlation strength between the 

discrete labels fi and v here used to describe the finite-dimensional phase space"? 

To answer this question, let us initially comprehend the role of J^{rj, £) present in 

the function E(jjl, v\\i' , v'): it was basically written as a sum of products of Jacobi 

^-functions evaluated at integer arguments, which plays, in such an aforementioned 

discrete phase space, the role reserved to the Gaussian functions in the continuous 

case (N —> 00); besides, y/4j + 2 represents its respective width and it assumes a 

constant value, in this situation, for a given j. Hence, any sub-Planck structured^ 

or even relevant correlations within this range are smeared in the smoothing process, 

which directly affects the description of correlation strength via Wehrl-type entropy 

functionals. Another important restriction associated with E[H; t] emerges from the 

mathematical property S[p(t)] = — Tr [p(i) mp(i]l < E[H;t], i.e., Eq. (1231) consists 

of an upper bound for the von Neumann entropy consequently, "if there are small 

distance fluctuations or if Jf p ([J,, v\ t) is concentrated on small regions of discrete 

phase space", then the percent error estimated S = S/E will be very bad. However, 

these apparent li mita tions can be circumvented by: (i) adopting the prescription of 

Manfredi and FebJ^ for quantum entropy based on Wigner functions in continuous 

phase-space, (ii) reformulating the function Jf(r], £) in order to modify its respective 

Fool \f\0\ 
width, and finally (iii) including a parallel study on von Neumann entropy^ and 
J7I 

quantum discord which permits us to increase our knowledge base on the intricate 
mechanisms of correlation strength related to the spin-squeezing and entanglement 
effects for any spin systems. In many ways the analysis presented in this work of the 
modified LMG model is complementary to that provided in Ref. [29J moreover, in 
what concerns the discrete phase-space approach and its implications for different 
physical systems, our results sound promising at a first glance. 

5. Conclusions 

In this non-trivial quantum mechanical scenario of correlations, entanglement and 
spin-squeezing effects, as well as their connections, any associated theoretical and/or 
experimental proposals for measures of the aforementioned effects shall necessarily 
be accomplished by exhaustive tests of confidence within a wide class of analogous 
physical systems. The mere comprehension of these effects by means of a 'specific 
theoretical/experimental measure applied to a particular physical system' does not 
change its respective status of proposal per se for standard measure: such transition 
demands 'time, patience, and efforts' to understand the subtle role of correlations 
and their intrinsic mechanisms in quantum mechanics. In this paper, we have used 
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a theoretical framework of finite-dimensional discrete phase space as an alternative 
approach for the study of the spin-squeezing and entanglement effects. The emphasis 
on covariances and its fundamental role in the investigation of spin-squeezing effects 
is not accidental: additional correlations related to the anticommutation relations of 
the angular-momentum generators are now included in the analysis of spin systems. 
This first effective gain does not imply in the match between the previous effects, 
since both the measures here adopted for describing entanglement exhibit 'slightly 
different functional relations' and also present 'speculative features', as expected. 

In what concerns the connection between the spin-squeezing and entanglement 
effects, both time-evolution operator and initial state of a given spin system play a 
fundamental role in such process: they describe the short-range and/or long-range 
correlations of the multipartite system under investigation, which could explain the 
similarity degree of the spin-squeezing and entanglement measures studied in this 
paper, as well as its 'almost perfect match'. In this sense, the modified LMG model 
fulfils a reasonable set of mathematical and physical prerequisites that lead us to 
corroborate, by means of numerical computations, the aforementioned link, for then 
establishing, subsequently, its inherent limitations. Besides, these results open new 
possibilities for future investigations in different physical systems which encompass 
underlying su(2) structures with distinct physical properties — in particular, those 
multipartite systems where contributions related to the long-range correlations are 
indeed necessary in the effective description of entanglement. 

Our particular phase-space approach, nevertheless, also reveals its drawbacks. 
To begin with, Eq. (|12l) and its constituting blocks coming from the numerical recipe 
described in subsection 3.2 do not satisfy the basic criterion 'easy-to-compute'. The 
large number of sums that appeared in the time evolution of discrete Wigner func- 
tion and mean values obligatorily implies in high computational and operational 
costs, whose complexities grow as N increases. Moreover, the modulo N extraction 
phase here adopted for the discrete variables still follows that mathematical pre- 
scription discussed in Ref. 1191 which represents a necessary operational cost inherent 
to the mod(7V)-invariant operator basis G(/x, v). However, such apparent limitations 
can be circumvented in this context by adopting the theoretical framework exposed 
in Ref. 24 for the mo d(iV)- invariant unitary operator basis A(/i, v) 1 and reformulat- 
ing the content associated with the time evolution of the discrete Wigner function. 
This procedure will allow a real computational gain in the numerical calculations. 

Now, let us briefly mention some possibilities for future research that stem from 
the present paper. As a first example, we recall from Ref. |7| those considerations on 
the concept of quantum discord and its important link with the subtle boundary 
between entanglement and classical correlations. Given that the Wehrl-type entropy 
functionals failed in the description of correlations related to the entanglement and 
spin-squeezing effects, it seems reasonable to introduce such a measure in this con- 
text since quantum discord can lead us to a more efficient analysis on quantum and 
classical correlations in multipartite physical systems. Another possible example of 
research consists in employing that finite-dimensional discrete phase-space frame- 
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work in order to corroborate the Lieb's conjecture for the spin coherent states^ 
via quantum dynamics of spin systems. 
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Appendix A. Variances, covariances, and uncertainty relation 

There are interesting results relating covariance functions and variances associated 
with certain summations of two non-commuting observables X and Y, namely, 

*x± Y = Vx + *V ± 2fxY and Y x±iY = fx - Y Y ± 211^. 

Note that such connections can be generalized for a large number of non-commuting 
observables. For instance, if one considers the set {X, Y, Z}, it is easy to show that 

r x+ Y+z - r x + r Y + r z + 2r XY + 2r xz + 2r YZ 

= ^X+Y + ^X+Z + ^Y+Z — ^JC — — ^Z 

yields the statistical balance equation 

*x+y+z + fx + "Vy + = Vx+y + Vx+z + Vy+z- (A.l) 

This particular 'conservation law' describes, if applied to the angular-momentum 
generators, how the statistical fluctuations are distributed among the observables 
{Jz, J y, J z } in the measurement process for a given initial quantum state. Table [2] 
shows some analytical expressions obtained for the variances and covariances related 
to the angular- momentum generators and spin coherent states, which exemplify Eq. 
(jA.lj) . Following, let us discuss how the unitary transformations modify the signal- 
to-noise ratio for an ideal experimental sit uat ion where the deleterious effects of a 
low-efficiency detection process are absent P^l 

The generator of unitary transformations T(£,u;) maintains, by definition, the 
linearity of the angular-momentum operators. Indeed, the new set 

{J a := Tt(£, w )J T(&w) : £ G C and w e ^\ a=x ^ z 

shows explicitly such a mathematical property through the relations 

Ja; = A\V&x + A^Jy + Ai 3 3 z 
J y = A2\i x + A22^ y + A23JZ 

Jz = A 31 3 x + A 32 J y + A 33 3 z , (A.2) 

whose coefficients {Ay} can be immediately determined from the results obtained in 
subsection 2.1. Since the commutation relation [J a , Jt] = ie a fc c J c keeps its invariable 
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Tabic 2. The explicit results for the variances and covariances — as shown on the 
table below — allow us, in principle, to illustrate both the Eqs. (|6jl and IjA.lj l, Besides, 
they also lead us to verify the uncertainty relation ¥j a fj b — (^J a J b ) > k | ([3a, Jj,] ) | 2 
for any a,b = x, y, z. Note that the saturation is reached in this case for 8 = and 7r, 
which correspond to the fiducial states \j, —j) and 

Variances and covariances associated with {3 x ,3 y ,3 z } for the spin coherent states 

n. =072) [l-cos 2 (v>)sin 2 (0)] 
^=072) [l-sin 2 (^)sin 2 (0)] 
Yj, =(i/2) sin 2 (0) 

=-074) 8111(299) sin 2 (0) 
*J,J, = 0/4)cos( V )sin(2e) 
ft y 3. = (j/4) sin(vp) sin(2e) 
*1,+J, = 0/2) {2 - sin 2 (e) [1 + sin(2^)]} 
*J«+J, = 0/2) [1 + sinV) sin 2 (9) + cos(^) sin(20)] 
*J B +J, = 0/2) [1 + cos 2 (<p) sin 2 (0) + sin(<p) sin(20)] 
*J*+J S +J, = 072) I 2 " sin(2 V ) sin 2 (0) + [cos(p) + sin((p)] sin(20)} 
Sec Rcf. 35 for certain mean values involving some quadratic forms of 3± and J z . 

form, some restrictions on the aforementioned coefficients should be established in 
this context, that is, 



A22A33 


— A23A32 


= A n , 


A13A32 


— A12A33 


= A 2 i, 


A12A23 


— A13A22 


= A31 


A23A31 


— A21A33 


= A 12 , 


A11A33 


- A 13 A 31 


= A 22 , 


Al 3 A 2 i 


-A n A 2 3 


= A32 


A21A32 


— A22A31 


= A13, 


^12^31 


- A n A 3 2 


= A23, 


A11A22 


-A 12 A 21 


= A33 



It is worth stressing that these general equations are valid for any unitary transfor- 
mations which preserve Eq. (|A.2[) ; moreover, after certain proper manipulations of 
such equations, this set can be reduced to 

A n A 12 + A21A22 + A 31 A 32 = 
A n A 13 + A21A23 + A 31 A 3a = 

A 12 A 13 + ^22-423 + A32A33 = 0. 

As mentioned in the main part of the text, the identity J = J 2 also preserves the 
original form of the total spin operator, which implies in the additional relations 

An + A 2 i + A 31 = 1, A\ 2 + A\ 2 + A 32 = 1, A 13 + A\ z + Al 3 = 1. 

Note that 'f-j + 7j = "Pj x + 'Pj v +'fj z = j represents an important by-product 

in this process since the sum of such variances remains invariant under the unitary 
transformation described by Eq. (|A.2[) . 

The evaluation of the signal-to-noise ratio (SNR) associated with collective and 
intrinsinc degrees of freedom is crucial in any experimental data analysis: indeed, it 
provides quantitative information on the measured signal and the noise inherent to 
the experiment under investigation. Thus, let us define SNR through the expression 
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£ft<j : = l(Ja)|/-\/^Ja' where each component of J presents an important role in the 
measurement process. However, if one considers experimental situations — rotations 
and/or time evolutions — that transform J a t— > J a , it seems natural to analyse the 
case d\ a i-> D\ a '■= KJn)|/\/"^I • ^ n * n ^ s sense ) the variances 

=A? 1 r J . + A? 2 ^ Jv + A? a r J , 

+2 (AuAi2t Us + ^ii ^i3^ j« + A 12 Ai 3 r Js jJ 

^ = ^21^ + ^2^ + ^23^* 

+ 2 (yl21^22^J x J v + ^2lA 23 rj x j z + ^22A 2 3^J y J 2 ) 

=^31^+^32^+43^ 

+ 2 (^ 3 lA32^J xJy + ^31^33^ J, + A 32 A 33 f Jyj!! ) 

exhibit an explicit dependence on all the previous variance and covariance functions, 
which modifies substantially the estimate of {£H a }- Besides, the uncertainty relation 

^ b -(Vljf>\\([Ja,J b })\ 2 (A-3) 
is also modified in this context, since the covariances 

y- y - iv = A 11 A 2l f Ja: + A 12 A 22 Y Jy + A 13 A 23 f Jz + {A n A 22 + A 12 A 21 )f J ^ y 
+(A n A 23 + A 13 A 21 )V^^ + (A 12 A 23 + A 13 A 22 )f iy ^ 
= A X1 A 31 1/^ + A 12 A 32 f Jy + A 13 A 33 Y^ + (A n A 32 + A 12 A 31 )%^ y 
+(A n A 33 + A 13 A 31 )f isc ^ + {A 12 A 33 + A 13 A 32 )f iy ^ 
Vj y j z = A 2l A 31 %^ + A 22 A 32 f Jy + A 23 A 33 Y 3z + (A 21 A 32 + A 22 A 3l )r^ 3y 
+(A 2 iA 33 + A 23 A 3 i)r^^ + (A 22 A 33 + A 23 A 32 )% yJz 

possess a similar dependence if one compares them with the previous case. Finally, 
it is worth stressing that such results are extremely relevant in the inves tigative 
process of squeezing and entanglement effects for spin systems JHHHIIHIII] 

For instance, Munoz and Klimov^ have recently introduced a particular set of 
discrete displacement generators upon a 2™ <E> 2™ discrete phase space, which allows 
us to establish a specific family of discrete spin coherent states for n-qubit systems 
with interesting mathematical and physical properties. Generated by application of 
the discrete displacement generators to a symmetric fiducial state, such states have 
isotropic fluctuations in a tangent plane whose geometric features are well-defined: 
in general, its direction does not coincide with that chosen for the mean value (J) . 
Besides, these reference states constitute the essential basic elements necessary for 
investigating the squeezing effects associated with the non-symmetric n-qubit states 
(resulting from the application of XOR gates) through the covariances Vj a j b - It is 
important to emphasize that the number of XOR gates applied to the discrete spin 
coherent states in order to minimize fluctuations in the original homogeneous plane 
strongly depends on the number of qubits. Since {^j a j,,} a ,b=x.y.z is also modified 
by the action of XOR gates, we believe that our results can be somehow useful for 
the detailed examination of such squeezing effects. 
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Appendix B. The Kitagawa-Ueda model 



Let us initiate our in vest igation on spin squeezing and entanglement through the 
Kitagawa-Ueda modeff^ — here described by the Hamiltonian H = x-J 2 — which 
describes a nonlinear interaction proportional to J 2 . The unitary transformations 
generated by the time-evolution operator U(r) = exp (— irJ 2 ) for r = %i represents 
a basic set of mathematical tools that allows us, in principle, to calculate certain 
important quantities necessary to the investigative process. Indeed, if one considers 
the unitary transformations 

J x (t) ee U^(t)J x U(t) = exp(ir) [J, cos(2tJ 2 ) - J y sin(2rJ,)] 
J v (t) ee U t (r)J y U(r) = exp(ir) [J x sin(2rJ 2 ) + 3 y cos(2rJ z )] 
J,(t) = Ut(r)J,U(r) =J Z ^ [J 2 , U(r)] = 0, (B.l) 

it is not so hard to yield the mean values (Heisenberg picture) 

(J x (t)) 9 , v = (0,<p\J x {t)\8,<p) =^ j - 1 jcos[cp-(2j-l)6 T }sm(e) 
(3 y (r))0, v ee (6,<p\J y (T)\d,<p) = */V- l j S m[(p- (2j - 1)6 T ] sm(d) 
(J z (r)) e>V3 = (6, tp\3 z (r)\9, <p) = -j cos(d) (B.2) 

i 

for the spin coherent states, where gf T = [cos 2 (r) + sin 2 (r) cos 2 (#)] 2 denotes the 
amplitude function, and 5 T — arctan [tan(r) cos(^)] the phase function. Note that 
such results lead us to conclude that unitary transformations involving nonlinear 
forms of the angular-momentum generators do not preserve, in general, the linear 
feature of such operators — for instance, see Eqs. (|B.1[) and (IA.2I) . This particularity 
is directly associated with the generators of the su(2) Lie algebra and does not apply 
to the squeeze operator S(() := exp [~ (C*a 2 — C a ^ 2 )] > since S(£) is responsible for 
unitary transformations in the Heisenberg- Weyl a lgebra which preserve the linearity 
of the boson annihilation and creation operatorsPZEH 

Another interesting point inherent to the time-evolution operator U(r) establi- 
shes the relation (3 2 (T))g^ — (J 2 (0))e i¥ , = + which preserves the mean value 
of the total spin operator in r = 0. In fact, this result can be interpreted as a direct 
consequence of the unitary transformations originated from the action of U(r) on 
the angular-momentum operators. Following, let us discuss some important points 
related to uncertainty relation, squeezing effects and their links with entanglement, 
considering the spin coherent states as initial states of the physical system. 



Initially, let us introduce the parameters S x := yj^/TZ^z and S y := fj y /TZ xyz 

i 

with TZ xyz := [('^j x j h ) 2 + j|(Jz)| 2 ] 2 (in this definition, lZ xyz ^ represents a 
condition sine qua non), which lead us to rewrite Eq. ([5]) as follows: S x S y > 1. 
In principle, this simplified form allows us to investigate the squeezing effects re- 
lated to the aforementioned model through exact expressions obtained in Table 
Ofor the time-dependent variance and covariance functions. Figure [TJa) shows 
the plots of S x (dot-dashed line) and S y (dashed line) versus r £ [0, 27r] for 
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Table 3. The time-dependent variance and covariance functions, as shown on the table below 
for the spin coherent states, constitute an important group of formal mathematical results related 
to the Kitagawa-Ueda model that permits us to investigate not only the Robertson-Schrodinger 
uncertainty principle, but also both the squeezing and entanglement effects associated with the 
collective angular-momentum operators. 

Time-dependent mean values, variances and covariance function for the Kitagawa-Ueda model 



<^MW=j/2 + 074)(2j-l) 
(3 2 y (r)) e , 9 =j/2 + (j/4)(2j-l) 



sin 2 (9) 
sin 2 (6») 



1 + ^2r~ 2 COS l 2 f - ( 2 i - 2 ) 5 2r] 

1 - cos [ 2 <P - (2J - 2 )<W] 

(J 2 z (r)) e , v = (J2(0)) e , v = j 2 cos 2 (0) + (j/2)sin 2 (0) 

y 3x (r) = j/2 + (j'/4)(2j - 1) [l + cos ^ ~ (2J - 2)*h-]] sin 2 (0) 

-j 2 ^T j ~ 2 cos 2 [tp - (2j - 1)S T ] sin 2 (6») 
(t) = j/2 + 0/4)(2j - 1) [l - sf^T 2 cos[2^ - (2j - 2)<5 2 r]] sin 2 (0) 

-j 2 ^r j ~ 2 sin 2 [ip - (2j - l)<S r ] sin 2 (6») 
r^(r) = '^(0) = (j/2)sin 2 (0) 

*J*J y to = 074)(2j - l)^ 2 2 T j - 2 sin[2 V - (2j - 2)8 2t ] sin 2 (0) 
~(j 2 /2)si/^~ 2 sin[2</> - (4j - 2)<5 T ] sin 2 (6») 

See Refs. llOl and H"2l for some measurement criteria involving such physical quantities. 



j = 2 and 8 = tp = j fixed. The hachured area exhibited in the picture descri- 
bes, in this case, that region where the squeezing effects occur for S x ^(t) < 1 
(but not both simultaneously). Such a particular evidence of squeezing effect 
gene rali zes, through an effective way, those results obtained by Kitagawa and 
UedJ^ in the absence of one-axis twisting mechanising] 
• How the squeezing and entanglement effects can be connected? S0rensen et aP^l 
have proposed an interesting experiment for the Bose-Einstein condensates (and 
attainable with present technology) which allows us, in principle, to answer this 
intriguing question through a fundamental effect in quantum mechanics: the 
many-particle entanglement. Since correlations among spins yield the squeezing 
effect, it is natural to build a bridge connecting both the entanglement and spin 
squeezing effects by means of a solid mathematical framework which allows us to 
establish certain reliable measurement criteria!^ In this sense, let us comment 
some few words on the entanglement criterion adopted in such an experiment: 
it is basically focussed on the separability criterion for the TV-particle density 
operator involving the variances and squared mean values associated with each 



g In fact, this mathematical procedure introduces artificially additional quantum correlations (by 
means of time-dependent unitary transformations that consist of rotations around the x-axis) for 
the angular-momentum operators, which allow to reduce the standard quantum noise related to 

the spin coherent states down to i V3 • It is important to emphasize that such an estimate does 
not consider the contributions originated from the covariance function. 
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Fig. 7. (a) represents the plots of S x (dot-dashed line) and S y (dashed line) versus r £ [0, 2tt] 
with j = 2 (AT = 4 spins) and = <p = fixed. Note that both the parameters 5a; (r) and S y (r) 
assume values less than one in different intervals of r (see hachured region), which implies in the 
squeezing effect for each case obeying the relation Sx(t)S v (t) > 1 (the saturation is reached for 
r = mn with m S N). (b) corresponds to the plots of £ x (dot-dashed line), £ y (dashed line), and 
£ z (solid line) versus r £ [0, 27r], for the same values of {j, 0,<p} adopted in the previous picture. 
The hachured area denotes, in this case, that region where the entanglement effect occurs. 



orthogonal component of the collective angular-momentum operators, that is, 



(J b ) 2 + (J c ) 2 



< 1 (a,b,c = x,y,z). 



(B.3) 



Thus, the parameters {£ a } characterize the atomic entanglement in this con- 
text, and quantum states with £ a < 1 are referred to as spin squeezed statesff^ 
For instance, Figure |7£b) shows the plots of £ x (dot-dashed line), £ v (dashed 
line), as well as £ z (solid line) versus r e [0, 2ir] for the same values of {j, 6, ip} 
adopted in the squeezing criterion — see Fig. EJa). Note that £ x and £ v ex- 
hibit values less than one, except £ z , which implies in entanglement effect (see 
hachured region) for the orthogonal components 3 X (r) and 3 y (r) , as expected, 
since [H, J z ] = 0. The similarities between both the figures and the evidence of 
the squeezing and entanglement effects for the same values of r justify, in such 
a case, the above connection per se. 



Appendix C. A case study for h — 

To what extent the anisotropy parameter 7 can be effectively employed in order to 
validate the entanglement criteria here adopted in the modified LMG model? Since 
the two-body term of the Hamiltonian H' = —hJ z —jj (J 2 , + 7J 2 ,) plays an essential 
role in such a case, we cannot refrain from studying it per se in a deeper perspective 
— in this context, we will assume h = hereafter. Note that the plethora of results 
exposed and discussed along the text for assigning the existence of entanglement 
in a multi-spin system can then be directly compared in the particular case of the 
two-body correlations associated with the model at hand. This important procedure 
will produce an effective range for 7 with N = 20 fixed, via extensive numerical 
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(a) (b) 




Time Time 



Fig. 8. Plots of (a,c) £ z (t) (solid line), E z (t) (dot-dashed line), S { z x > (t) (dashed line) and (b,d) 
3. (t) versus t G [0, 50] for (a,b) 7 = 0.5 and (c,d) 7 = 0.948 (critical value), with N = 20 fixed 
and h = (absence of transverse magnetic field). The transition between the different values of 
7, here depicted by pictures (a) and (c), illustrates how the squeezing and entanglement criteria 
(as well as the covariance function) are affected by the anysotropy parameter 7. Such a numerical 
evidence suggests the existence of a validity domain for the aforementioned criteria, where the 
'almost perfect match' between squeezing and entanglement effects should be preserved for any 
t > 0. 



computations of Eqs. (|B.3|) and (|2"5)l , whose physical implications permit, within 
other features, to answer (in part, at least) the aforementioned question. 

Now, we discuss some essential points inherent to the numerical results obtained 
through the discrete Wigner function (fT3|) evaluated at equally spaced time intervals 
(i.e., 0.05 time units) and I7I < 1. The first one confirms the link between squeezing 
and entanglement effects verified in section 4 for values up to 7 = 0.5; however, if 
one considers the interval 0.5 < 7 < 0.9, small differences concerning the parameters 
{Sa(t),£ a (t), E a (t)} for a — y,z begin to appear, this fact being considered by us 
as a first numerical evidence about the 'limitation' of the particular entanglement 
measures £ a (t) and E a {t) here studied. The second essential point establishes that, 
for 7 0.948 (critical value), £ a (t) does not predict entanglement effects anymore, 
while E a (t) still does — in this case, it is worth stressing that S a x \t) also indicates 
the presence of the squeezing effect for determined values of time. 

With respect to the proposals of squeezing and entanglement measures discussed 
in this work, let us draw attention to the presence of the covariance function "V3 j z (t) 
in the expression for Sa (t) (see table [1] and appendix A) which depends explicitly 
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on the mean value (|18[) . Our numerical calculations allow to show, in particular, that 
"fj v j,{t) is responsible for the occurrence of squeezing effect in the entire domain of 
7, since it carries essential additional information, when we take into account the 
RS inequality (instead of the Heisenberg one) , associated with the anticommutation 
relation of the angular-momentum operators J y and J z . 

Therefore, our results reinforce the real necessity of a deeper discussion on the 
validity domain related to the entanglement criteria here used to characterize the 
important link between squeezing and entanglement effects!^ Thus, it seems that 
we cannot assure beforehand the viability of such criteria without also especifying 
their respective validity domains in order to be employed with a consistent degree 
of confidence. 
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